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Abstract 

We study two-bodies gravitational problem where the mass of one of the 
bodies varies and suffers a damping-anti damping effect due to star wind 
during its motion. A constant of motion, a Lagrangian and a Hamiltonian 
are given for the radial motion of the system, and the period of the body is 
studied using the constant of motion of the system. An application to the 
comet motion is given, using the comet Halley as an example. 
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1 Introduction 



There is not doubt that mass variable systems have been relevant since the founda- 
tion of the classical mechanics and modern physics too [0] which have been known as 
Gylden-Meshcherskii problems [1]. Among these type of systems one could mention: 
the motion of rockets [2], the kinetic theory of dusty plasma [3], propagation of elec- 
tromagnetic waves in a dispersive-nonlinear media [4], neutrinos mass oscillations 
[5], black holes formation [6], and comets interacting with solar wind [7]. This last 
system belong to the so called "gravitational two-bodies problem" which is one of 
the most studied and well known system in classical mechanics [8]. In this type of 
system, one assumes normally that the masses of the bodies are fixed and unchanged 
during the dynamical motion. However, when one is dealing with comets, beside to 
consider its mass variation due to the interaction with the solar wind, one would 
like to have an estimation of the the effect of the solar wind pressure on the comet 
motion. This pressure may produces a dissipative-antidissipative effect on its mo- 
tion. The dissipation effect must be felt by the comet when this one is approaching 
to the sun (or star), and the antidissipation effect must be felt by the comet when 
this one is moving away from the sun. 

In previous paper [14]. a study was made of the two-bodies gravitational prob- 
lem with mass variation in one of them, where we were interested in the difference 
of the trajectories in the spaces (x,v) and (x,p). In this paper, we study the two- 
bodies gravitational problem taking into consideration the mass variation of one 
of them and its damping-anti damping effect due to the solar wind. The mass of 
the other body is assumed big and fixed , and the reference system of motion is 
chosen just in this body. In addition, we will assume that the mass lost is expelled 
from the body radially to its motion. Doing this, the three-dimensional two-bodies 
problem is reduced to a one-dimensional problem. Then, a constant of motion, the 
Lagrangian, and the Hamiltonian are deduced for this one- dimensional problem, 
where a radial dissipative-antidissipative force proportional to the velocity square 
is chosen. A model for the mass variation is given, and the damping-anti damping 
effect is studied on the period of the trajectories, the trajectories themselves, and 
the aphelion distance of a comet. We use the parameters associated to comet Halley 
to illustrate the application of our results. 
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2 Equations of Motion. 



Newton's equations of motion for two bodies interacting gravitationally, seen from 
arbitrary inertial reference system, and with radial dissipative-antidissipative force 
acting in one of them are given by 

d ( dx\\ Gm\m 2 , . M . 

m i-rr = —, n(ri-r 2 ) (la) 



dt \ dt J Iri — r 2 | 3 



and 



d ( dr 2 \ Gmim 2 , . 7 

"77 rn 2 — = --, p7(r 2 - r x ) - 



dt \ dt J |r 2 — ri I 3 Iri — r 2 | 



d|ri - r 2 
dt 



(r 2 - n) , (lb) 



where mi and r«2 are the masses of the two bodies, ri = (x\, yi,Zi) and 
i"2 = (x2,V2, 22) are their vectors positions from the reference system, G is the 
gravitational constant (G = 6.67 x 10~ n m 3 / Kg s 2 ), 7 is the nonnegative constant 
parameter of the dissipative-antidissipative force, and 



|ri - r 2 | = |r 2 - r x | = v 7 i x 2 ~ ^i) 2 + (2/2 - Vi) 2 + fa ~ 
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is the Euclidean distance between the two bodies. Note that if 7 > and 
d\ri — r 2 |/dt > one has dissipation since the force acts against the motion of the 
body, and for d\ri — r 2 |/dt < one has anti-dissipation since the force pushes the 
body. If 7 < this scheme is reversed and corresponds to our actual situation with 
the comet mass lost. 



It will be assumed the mass mi of the first body is constant and that the mass 
ni2 of the second body varies. Now, It is clear that the usual relative, r, and center 
of mass, R, coordinates defined as r = r 2 — ri and R = (miri+rr^i^) /(mi+7712) are 
not so good to describe the dynamics of this system. However, one can consider the 
case for ra\ ^> m 2 (which is the case star-comet), and consider to put our reference 
system just on the first body (r x = 0). In this case, Eq. (la) and Eq. (lb) are 
reduced to the equation 



m 2 



~dtf 



Gm\m 2 



r — m 2 r — 7 



dr 
dt 



(2) 



where one has made the definition r = r 2 = (x,y,z), r is its magnitude, 



or + y 



2 + z 2 , and r = r/r is the unitary radial vector. Using spherical coordinates 
x = r sin 9 cos ip , y = r sin 9 sin cp , 2; = rcos6*, (3) 
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one obtains the following coupled equations 

m 2 [r — ru — rip sin 9) — m 2 r — , (4) 

m 2 (2r9 + r9 — rip 2 sin 9 cos 9) = —rh 2 r9 , (5) 

and 

m 2 (2fip sin 9 + r(p sin 9 + 2r0^ cos 9) = —rfi2rip sin 9 . (6) 
Taking = as solution of this last equation, the resulting equations are 



m 2 (r — r9 ) = m^r — , (7) 

and 

m 2 (2r9 + r9) + m 2 r9 = . (8) 

From this last expression, one gets the following constant of motion (usual angular 
momentum of the system) 

k = m 2 r 2 9 , (9) 

and with this constant of motion substituted in Eq. 7, one obtains the following 
one-dimensional equation of motion for the radial part 

d 2 r Gnu m 2 fdr\_ + I 2 ^ 



dt 2 r 2 m 2 \dt J m 2 m 2 r 3 

Now, let us assume that m 2 is a function of the distance between the first and the 
second body, m 2 = m 2 (r). Therefore, it follows that 

m 2 = m' 2 r , (11) 

where m' 2 is defined as m' 2 = dm 2 /dr. Thus, Eq. (10) is written as 

d 2 r Gmi l 2 e m' 2 + 7 f dr\ 2 

' (12) 



dt 2 r 2 m 2 r 3 m 2 \dt 

which, in turns, can be written as the following autonomous dynamical system 
dr dv Gm l l 2 e m' 2 + 7 2 

Tt =v; Tt=-— + ^-^ 2 - v ■ (13) 

Note from this equation that m' 2 is always a non-positive function of r since it 
represents the mass lost rate. On the other hand, 7 is a negative parameter in our 
case. 
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3 Constant of Motion, Lagrangian and Hamilto- 
nian 



A constant of motion for the dynamical system (13) is a function K = K(r, v) which 
satisfies the partial differential equation [9] 



dK 
dr 



-Gmi Q m' 2 + 7 



+ -v 



r 2 m 2 r 3 vri2 



£-■>■ 



The general solution of this equation is given by [10] 

K(x,v) = F(c(r,v)) , (15a) 

where F is an arbitrary function of the characteristic curve c(r, v) which has the 
following expression 



c{r, v) = m 2 (r)e 



2 7 A(r V + f f ^Gm, _ 2% \ m 2 {r)e 2,X(r) dr ^ 

J \ r z m^r 6 J 



and the function A(r) has been defined as 

/dv 
m 2 (r) 

During a cycle of oscillation, the function rri2(r) can be different for the comet 
approaching the sun and for the comet moving away from the sun. Let us denote 
m 2 +(r) for the first case and wi2_(r) for the second case. Therefore, one has two cases 
to consider in Eqs. (15a), (15b) and (15c) which will denoted by (±). Now, if 
denotes the mass at aphelium (+) or perielium (-) of the comet, F(c) = c^j1m° 2± 
represents the functionality in Eq. (15a) such that for m 2 constant and 7 equal 
zero, this constant of motion is the usual gravitational energy. Thus, the constant 
of motion can be chosen as = c(r, v)/2m\±-, that is, 

^ ± = ^ e27A±(rV + ^/ (r) ' (16a) 
where the effective potential V e ff has been defined as 

eff[r ' m° 2± J r 2 m° 2± J r 3 [ ' 
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This effective potential has an extreme at the point r* defined by the relation 

/ 2 



r*m 2 (r* 



;i7) 



which is independent on the parameter 7 and depends on the behavior of 7712(7") 
This extreme point is a minimum of the effective potential since one has 



d 2 V ± 
dr 2 



> . 



Using the known expression [11-13] for the Lagrangian in terms of the constant of 
motion, 

K(r, v) dv 



L(r, v) = v 



(19) 



the Lagrangian, generalized linear momentum and the Hamiltonian are given by 

L± = (20) 



HMl e ^±(r) v 2 _ y± (r) 

2m° 2± effK ' 



V 



m 



2± 



and 



± _ m° 2± p 2 27A±(r) ± 



(21) 



(22) 



The trajectories in the space (x,v) are determined by the constant of motion (16a). 
Given the initial condition (r ,v ), the constant of motion has the specific value 



m l±( r o) 2 7 A ± (r ) 2 , y± / N 



(23) 



and the trajectory in the space (r, v) is given by 



/^^± e -7A ± (.) 



m| ± (r) 



^ - V5,(r) 



nl/2 



(24) 



Note that one needs to specify o also to determine Eq. (9). In addition, one 
normally wants to know the trajectory in the real space, that is, the acknowledgment 
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of r = r(6). Since one has that v = dr/dt = (dr/d8)6 and Eqs. (9) and (24), it 
follows that 



0(r) = o + 



m 2± 



(r)e 7A±w dr 



V 2m °2± Jr r 2 JK± -V e %(r] 



The half-time period (going from aphelion to perihelion 
be deduced from Eq. (24) as 



1 l/2 



m2±(r)e rX± ^ dr 



(25) 

or backward (-)) can 
(26) 



where r\ and r 2 are the two return points resulting from the solution of the following 
equation 



1,2. 



(27) 



On the other hand, the trajectory in the space (r, p) is determine by the Hamiltonian 
(21), and given the same initial conditions, the initial p Q and Hf are obtained from 
Eqs. (21) and (22). Thus, this trajectory is given by 



P=±A 



2mL(r) g7A±(r) 



■in 



2± 



■nl/2 



(28) 



It is clear just by looking the expressions (24) and (28) that the trajectories in the 
spaces (r,v) and (r,p) must be different due to complicated relation (21) between v 
and p (see reference [14]). 



4 Mass- Variable Model and Results 

As a possible application, consider that a comet looses material as a result of the 
interaction with star wind in the following way (for one cycle of oscillation) 



m 2-( r 2(i-i)) ( 1 — e ar J incoming {-V) v < 

m 2 ±(r) = \ 

m 2+ (r 2 j_ 1 ) — b^l — e - "^ - ^- 1 ^ outgoing^—) v > 



(29) 

where the parameters 6 > and a > can be chosen to math the mass loss rate in 
the incoming and outgoing cases. The index "i" represent the ith-semi-cycle, being 
r 2(i-i) and r 2 i-i the aphelion(r a ) and perihelion (r p ) points (r D is given by the initial 
conditions, and one has that m 2 _(r ) = m ). For this case, the functions A + (r) and 
A_(r) are given by 

A+(r) = — 



and 



A_(r) 



-1 



a(b — m p ) 



ar 



+ In (m p - 6(1 - e~ a{r - rp) )) 



(30a) 



(306) 



where we have defined m a = m 2 (r a ) and m p = m 2 (r p ). Using the Taylor expansion, 
one gets 



e 2-yX+(r) _ e 2-yr/m a 



^ e _ Qr+ l_2 7 ( 2 7 



am a 



2 am a \am a 



1 ) e~ 2ar + 



(31a) 



and 



3 2 7 A-(r) 



2-fr 
g (b-m p ) 



27 



(m p - 6) a(m f- 
27 



1 + 



27 



21 



-o(r— r p ) 



a(m p — 6) m p — b 



- 1 



2 a(m p — 6) \a(m p — 6) / (m p — 6) 



-2a(r—r p ) 



The effective potential for the incoming comet can be written as 



(316) 



Ve + ff(r) 



Gmim a l 2 e 1 



2m a r 2 



(32) 



and for the outgoing comet as 



Gm x m a ll 1 



2m a r 2 



27r 
g (mp-b) 



- + W 2 (7,a,r), (33) 
(m p - 6) "('"p- 6 ' 



where W\ and are given in the appendix. 
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We will use the data corresponding to the sun mass (1.9891 x 10 30 i^g) and the 
Halley comet [15-17], 

m c « 2.3 x 10 u Kg, r p « 0.6 au, r a « 35 aw, Z fl « 10.83 x 10 29 ^-m 2 /s, (34) 

with a mass lost of about <5m « 2.8 x lO 11 ^ per cycle of oscillation. Although, 
the behavior of Halley comet seem to be chaotic [18], but we will neglect this fine 
detail here. Now, the parameters a and "6" appearing on the mass lost model, Eq. 
(29), are determined by the chosen mass lost of the comet during the approaching to 
the sun and during the moving away from the sun (we have assumed the same mass 
lost in each half of the cycle of oscillation of the comet around the sun). Using Eqs. 
(32) and (33), Eq. (24), the trajectories can be calculated in the spaces (r, v) . Fig. 
1 shows these trajectories using 5m = 2 x 10 w Kg (or 5m jm = 0.0087%) for 7 = 
and (continuos line), and for 7 = —3 Kg/m (dashed line), starting both cases from 
the same aphelion distance. As one can see on the minimum, dissipation causes to 
reduce a little bit the velocity of the comet , and the antidissipation increases the 
comet velocity, reaching a further away aphelion point. Also, when only mass lost is 
considered (7 = 0) the comet returns to aphelion point a little further away from the 
initial one during the cycle of oscillation. Something related with this effect is the 
change of period as a function of mass lost (7 = 0). This can be see on Fig. 2, where 
the period is calculated starting always from the same aphelion point (r a ). Note 
that with a mass lost of the order 2.8 x lO 11 ^^ (Halley comet), which correspond to 
Sm/m = .12%, the comet is well within 75 years period. The variation of the ratio of 
the change of aphelion distance as a function of mass lost (7 = 0) is shown on Fig.3. 
On Fig. 4, the mass lost rate is kept fixed to 5m/ m = 0.0087%, and the variation of 
the period of the comet is calculated as a function of the dissipative-antidissipative 
parameter 7 < (using I7I for convenience). As one can see, antidissipation always 
wins to dissipation, bringing about the increasing of the period as a function of 
this parameter. The reason seems to be that the antidissipation acts on the comet 
when this ones is lighter than when dissipation was acting (dissipation acts when 
the comet approaches to the sun, meanwhile antidissipation acts when the comet 
goes away from the sun). Since the period of Halley comets has not changed much 
during many turns, we can assume that the parameter 7 must vary in the interval 
(—0.01, 0]Kg/m. Finally, Fig. 5 shows the variation, during a cycle of oscillation, of 
the ratio of the new aphelion (r' a ) to old aphelion ( function of the parameter 

7- 
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5 Conclusions and comments 



The Lagrangian, Hamiltonian and a constant of motion of the gravitational attrac- 
tion of two bodies were given when one of the bodies has variable mass and the 
dissipative-antidissipative effect of the solar wind is considered. By choosing the 
reference system in the massive body, the system of equations is reduce to 1-D 
problem. Then, the constant of motion, Lagrangian and Hamiltonian were obtained 
consistently. A model for comet-mass-variation was given, and with this model, a 
study was made of the variation of the period of one cycle of oscillation of the comet 
when there are mass variation and dissipation-antidissipation. When mass variation 
is only considered, the comet trajectory is moving away from the sun, the mass lost 
is reduced as the comet is farther away (according to our model), and the period of 
oscillations becomes bigger. When dissipation-antidissipation is added, this former 
effect becomes higher as the parameter 7 becomes higher. 

It is important to mention that if instead of loosing mass the body would had 
winning mass, the period of oscillation of the system would decrease. One can 
imagine, for example, a binary stars system where one of the star is winning mass 
from the interstellar space or from the other star companion. So, due to this winning 
mass, the period of the star would decrease depending on how much mass the star 
is absorbing. 
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6 Appendix 

Expression for W\ and W2' 



Gm 2 \ v(v — l)e(~ 4:+p ' >ar 
W l = — ^{ J - + apEi{apr) - 2ap(p - l)#i((-4 +p)ar) 



in 



2+ 



2r 



ap 2 (p — 1) 



Ei((-A + p)ar) 



P<J) X) r ^3a(l -p)rE i ({p-3)ar)} 



p(p + 3 ) 



2 

p + 2 



3 (p-2)ar 



+ ar(p-2)£ i ((p-2)ar) 



2m2 + r 2 



a 2 p(p- l)r c(p _ 2)ar 



,(p-2)ar _ p e (v- l ) ar 



ap(p — l)e 



(p— 2)ctr 



ap(p — 1) 



2 

+ pare (p " 1)ar - p 2 are (p - 1)ar - p 2 a 2 r 2 E l (par) 



a 2 (p — 2) 2 p(p — l)r 5 



Ei((p — 2)ar) + pa 2 r 2 Ei((p — l)ar) 



-2a 2 p 2 r 2 Ei((p - l)ar) + p 3 a 2 r 2 Ei((p - l)ar) 



(Al) 

where m a is the mass of the body at the aphelion, and we have made the definitions 

2 7 



P 



ftm„ 



(A2) 



and the function is the exponential integral 

EAz) = 



00 e - t 



— dt 
t 



(A3) 
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Wo 



Gm% \ e^ c 



TTln 



i + J^-V^ + -JL 

2{m p + aq) m p + aq 



D gar 



+qaE l {qar)- q{q f, 



(m p + aq) 2 r 



ge 



gar 



(m p + ag)r 



,(q-3)ar 



a(g - 3)r^i((g - 3)ar)] - 2aEi((q - 2) or) 



+aqE i ((q - 2)ar) - / ^((g - 2)ar) 

+^ - 2)ar) ■ E % ((g - 2)ar) 



2(m p + ag) 2 
2g 2 ae^ r 



m„ + aq 



Ei((q - 2)ar) + - [ e ^~ x > r - (g - l)ar^((g - l)ar)] 



(m p + aq)r 
np qar 

f^- r [e (g - 1)ar - (g - l)ar^((g - l)ar)] 

[m p + aq)r 

I 2 



qae 



qar 



2m 2 + {m p + aq) c - 



+ q 2 a 2 Ei[qar) 



+ m , ' ^ n [-1 + 2ar - qar + (2 - g) VrV 2 -"^ ((g - 2)ar)] 



2{m p + aq) 2 r 
q e ( 2 Q~ 1 ) ar 



(m p + aq)r 2 



[-1 + ar + qar + (g - \fa 2 r 2 e^- q)ar E i ((g - l)ar)] 



(AA) 

where m p is the mass of the body at the perihelion, and we have made the definition 

2l 



q 



a{m p — b) 



(A5) 
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Figure 1: Trajectories in the (r,v) space with 5m/m = 0.009. 




o » n q u^i czj 



Figure 2: Period of the comet as a function of the mass lost ratio. 
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Figure 3: Ratio of aphelion distance change as a function of the mass lost rate. 




Figure 4: Period of the comet as a function of the parameter 7. 
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Figure 5: Ratio of the aphelion increasing as a function of the parameter 7. 
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